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LONG-RANGE INTERACTION POTENTIALS FOR OXYGEN 


By Willard E. Meador 
Langley Research Center 

SUMMARY 

Electron configurations and valence bond interaction energies are derived for the 
18 molecular combinations of two ground-state oxygen atoms. The method differs from 
previous valence bond calculations in three important respects: (1) correlations between 
molecular and atomic states are included, (2) proper relations between the various multi- 
electron bonds and antibonds are employed, and (3) configuration interactions between 
molecular states of the same symmetry are determined. Many of the long-range inter- 
action potentials are affected substantially by these modifications. 

INTRODUCTION 

Long-range interaction forces between atoms are significant input data for calcu- 
lations of the nonequilibrium (transport) properties of collision-dominated gases. Often 
used for this purpose is the approximation of perfect pairing (ref. 1), in which molecular 
energies are expressed as linear combinations of bond and antibond energies based upon 
symmetry and spin pairing considerations. Although the absolute results for individual 
collisions may be inaccurate, the method is thought to establish fairly reliable and useful 
relations between the various potentials operating in a gas. (See ref. 2 for a further dis- 
cussion.) Hence, if a small number of such potentials are available from ab initio 
quantum mechanical calculations or experimental data, the others will follow from the 
simple perfect pairing expressions. 

A principal objective of the present research is to improve the existing relations 
(ref. 3) between the 18 long-range interaction potentials experienced by colliding ground- 
state oxygen atoms. The present paper is thus a continuation of reference 4, which 
reviewed and analyzed the conventional perfect pairing theory and applied it, with modifi- 
cations, to the relatively simple example of molecular nitrogen. Also found in reference 4 
is the assertion that previous treatments (ref. 3) of grazing oxygen collisions are ques- 
tionable for three reasons: (1) the electron configurations of molecules and atoms were 
not correlated to prevent dissociation products from being mixtures of several atomic 
states, (2) the strengths of the various multi electron bonds and antibonds were not prop- 
erly related according to overlap and other considerations, and (3) configuration 



interactions between molecular states of the same symmetry were not included. These 
deficiencies are corrected in the current research. 

The procedure is outlined as follows: the determination of the nine wave functions 
pertinent to oxygen 3p atoms, the deduction of molecular electron configurations from 
correlation studies of the atomic states, and the calculation of interaction potentials by 
valence bond methods applied to the electron configurations. Interactions between molec- 
ular states of the same symmetry are imposed automatically by the correlation analysis, 
as are the magnitudes and directions of the corresponding energy shifts. All 18 inter- 
atomic potentials are found to differ from the expressions in reference 3, the trend defi- 
nitely being toward simpler representations, more convenient relations between poten- 
tials, and increased symmetry with respect to each other and the dissociation asymptote. 

Perhaps the most pronounced effect of the correlation analysis is the flattening of 
nine of the 18 potential curves in their outer regions. Hence, the contributions of these 
states to the collision integrals and cross sections of kinetic theory are somewhat dimin- 
ished from previous calculations, which, therefore, appear to underestimate the transport 
coefficients of atomic oxygen at high temperatures. Some confirmation of this flatness 
is obtained for the ground state of molecular oxygen by extrapolating the corresponding 
Rydberg-Klein-Rees semiempirical potential curve to large internuclear separations. 

The present results for this particular example are more easily connected to the extrap- 
olated semiempirical curve than are the results of reference 3. 

Finally, extrapolations of selected long-range electron structures and potential 
energies to intermediate internuclear separations successfully predict the ordering of 
the six most stable molecular states and suggest a number of potential humps and curve 
crossings in the process. 


SYMBOLS 


a,b atoms or molecular states 

E equilibrium energy 

g± molecular orbital eigenstates of electron angular momenta 

H a b interaction matrix element between molecular states a and b 

J a ,J v two-electron exchange energies of a- or 7 r-type symmetry; sometimes 

denoted J 2 a and J2 n , or simply J when speaking figuratively 
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JlcpJl it one-electron exchange energies of o- or 77 -type symmetry 

Ml magnetic quantum number for total electronic orbital angular momentum 

(atom) 

M s magnetic quantum number for total electron spin (atom) 

P x >Py>Pz atomic p functions or orbitals 
Q collision cross section 

Ql one-electron Coulomb energy 

u* molecular orbital eigenstates of electron angular momenta 

x,y,z Cartesian coordinates, also abbreviated notation for atomic orbitals p x , Py, 

Pz 


electron spin functions 

A magnetic quantum number for total electronic orbital angular momentum 

(molecule) 

<p atom-atom interaction potential 

<p' eigenvalue solution of secular energy equation 


<P 3>4>4 


potentials for one-, two-, three-, and four-electron bonds and antibonds 


<£ a ,</>b interaction potentials for molecular states a and b of same symmetry 
(no configuration interaction) 

^a^b same as <£ a and cp ^ but with configuration interaction 

Subscripts: 


Tt gerade M or even symmetry with respect to geometric center of oxygen 
molecule 


3 



u 


"ungerade" or odd symmetry with respect to geometric center of oxygen 
molecule 


Special notations: 

3 P, 1 D states of atomic oxygen 


^Ag, A 1 ITg, . . . standard molecular states 


^g’^U’^g’l 

^u’VV 


1 , 1 , 0 


(Ho|i) 

1,2, 3, 4 


+>- 


standard molecular orbital notations; 7 T-orbitals sometimes appear without 
superscripts 

atomic orbital combinations x + iy, x - iy, z, respectively 
typical atomic state function defined by equation (1) 
angular momentum-spin state for atom 
sets of coordinates for elections 1, 2, 3, 4 

superscripts meaning even or odd symmetry with respect to plane containing 
internuclear axis (z) 

ATOMIC 3 P FUNCTIONS 


The p-electron configurations of the nine angular momentum-spin states 
for oxygen 3p atoms are fundamental to the construction, by correlation analyses, of 
the long-range electron configurations of the 18 molecular states dissociating into nor- 
mal products. A typical example is the (1,1) state corresponding to the atomic wave 


function 


(i i o|i) = (4.y 



KiMD 

1(1)0<1) 

0(1)q(1) 

1(1)13(1) 

(4.')-V2 

l(2)o(2) 

1(2)o(2) 

0(2) a^2) 

1(2)13(2) 

l(3)o<3) 

1(3)o<3) 

0(3)o<3) 

1(3)13(3) 


l(4)of(4) 

1(4)o(4) 

0(4)q<4) 

1(4)18(4) 


(i) 


where 1, 1, and 

quantum numbers 
arguments 1, 2, 

(lToJi). 


0 are one-electron functions pertaining to the individual magnetic 
1, -1, and 0, respectively, a and /3 are spin functions, and the 
3, and 4 refer to electrons. Equation (1) defines the notation 
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Two of the remaining eight ^Ml,M s J states are obtained immediately from equa- 
tion (1) by the application of the familiar ladder or step-down operator in angular momen- 
tum space. They are 

(0,1) : (lTOlo) ( 2 ) 

and 

(-1,1) : (l T 0|l) (3) 

However, configurations (1), (2), and (3) are not directly subject to the simple geo- 
metric interpretation required for correlation studies and the application of the approxi- 
mation of perfect pairing. A better representation for such purposes is obtained by 
writing 1 = x + iy, T = x - iy, and 0 = z, where x, y, and z stand for p x , p y , and 
p z atomic orbitals, respectively, and by expanding the wave functions in these terms. 

For example, 

|^1 1 0| lj = |(x+iy x-iy z | x+iy) 

= (x y z | x+iyj 

= (xy z|x) + i(x y z|y) (4) 

Similar expansions apply to all wave functions to yield the following complete set of 


oxygen 3p electron configurations: 

(1,1) : (x y z|x) + i(x y z|y) (5) 

(1,0) : (x y|x z) - (x z|x y) + i(x y|y zj - i(y z|x y) (6) 

(1,-1) : (x|x y zj + i^y |x y z) (7) 

(0,1) : (xyz|z) (8) 

(0,0) : (x z|y z) - ^y z|x zj (9) 

(0,-1) : (z|xyz) (10) 

(-1,1) : (x y z|x) - i(x y z|y) (11) 
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(-1,0) : (x y|x z) - ^x z|x y) - i(x y|y z) + i(y z|x y) (12) 

and 

(-1,-1) : (x|x y z) - i(y|x y z) (13) 

Expressions (6), (7), (9), (10), (12), and (13) are deduced from expressions (5), (8), and 
(11) by using the step-down operator in spin space. 

MOLECULAR ELECTRON CONFIGURATIONS 


Long-range molecular electron configurations are found by combining atomic states 
(5) to (13) in every possible way corresponding to two colliding oxygen 3p atoms a 
and b. If the z-axis coincides with that of the molecule, a p z electron on atom a 
combines with another on atom b to form either a (a) configuration for paired spins 
or a ( CT gj( a u) configuration for parallel spins; similarly, three and four p z electrons 

yield and (cr) 4 configurations, respectively, and ir*, ttJ, tt~, and -n~ 

orbitals are occupied by p x and p y electrons. The g and u subscripts standing 
for even and odd symmetries with respect to the molecular center are purposely absent 
in (cr)2 and (a ) 4 to indicate the valence-bond character of these configurations near 
the dissociation limit (ref. 4). A further discussion is given in appendix A. 

A simple example of the complete correlation procedure is the following develop- 
ment for the C >2 (^&gj state: 



(l,l) a + (l,l) b 


[(x y z|x) a + (x y z|x) b ] +[(x y z|y) a + (x y z|y) b ] 
+ [(x y z|x) a + (x y z|y) b ] + [(x y z|y) a + (x y z|x) b ] 


2“ 1 ( ff g)(^)[( ir+)4 W( ,r g) + (<)(4) (7r_)4 

/ _i_\ 2 / +\/ ~\/ — \2 / -A/ -f-\2 / -^2/ 

(^uy \ S/\ s) \J^m J \g) 


(14) 


Except for the last part of this expression, in which subconfigurations are defined from 
the previous part and each has a one-fourth probability of being measured, the positive 
signs within the various brackets and also in (l,l) a + (l,l) b denote schematically the 
interaction of atoms a and b. Finally, since the first two subconfigurations in the last 
part of expression (14) establish g-type symmetry for the molecule as a whole, the last 
two subconfigurations are chosen accordingly from several possibilities. 


6 



Although comparisons between expressions (5) and (11) show the atom-atom combi- 
nations in the middle part of expression (14) to be equally appropriate for (l,l) a + (-l,l) b 
collisions, the situation is complicated by the fact that £ molecular states of opposite 
symmetries with respect to a plane containing the z-axis are usually nondegenerate and 
must be so distinguished by + or - superscripts (ref. 5). Hence, the following results 

are obtained by arbitrarily suppressing either the pair [(x y z|x) a + (x y z|x) b ] and 

[(xy z|y) a + (xy z|y) b ] or the pair [(x y z|x) a + (x y z J y) J and [(x y z|y) a + (x y z|x) b ] 

in each collision: 


a 5 S+ : 2 " 1/2 ^g)f u )[( 7r+ ) 4 (%)( 7r i) + 


(15) 


5 S " : 2" 1 / 2 


(°g) ( G u) (<) ( ir g)( 7 r u)( ir g) + ( Wu ) (^g)' (*u) (^s) 




and 


5s u : 2 ^ ( a g)( a u) (0 (^(%) (^g) + (^u)^!) (^(^g) 


(16) 


(17) 


The validity of the arbitrary suppression of pairs of interaction terms in the devel- 
opment of expressions (15) to (17) is not obvious and must be investigated thoroughly. 
Three observations are pertinent to this objective: 

(1) Molecular A states cannot be written as linear combinations of £ states. 

Expressions (14) to (16) contradict this statement unless the a^Z)+ and states 

o o 

either do not occur or do not have the assigned electron configurations. The second 
alternative requires configuration interactions with additional states of the same 
symmetry. 

(2) The arbitrary suppression of pairs of interaction terms would not be required 

for the a ^g and states if the corresponding long-range energies were identical. 

Linear combinations of 5s states with different symmetries and equal energies are 
permitted by quantum theory if the resulting electron configurations are also ^ 2 . Such 
linear combinations would restore each pair of interaction terms to equivalent roles of 
prominence in the final electron descriptions. 

(3) The combination of (0,l) a with (0,1)^ to yield 

*> 5s g = < t ’> 4 ( 0 (’ g)(^)(’i) < 18 > 

is the only 62 possibility not yet considered. 
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Observations (1) and (3) imply that the 5 2 " state cannot be formed from ground- 
state oxygen atoms, nor can the oxygen molecule dissociate into normal atoms. 

This conclusion agrees with the more complex analysis based on group theory (ref. 5, 
p. 321) and, thus, supports the present correlation procedure. In addition, the a 5s+ 
state of expression (15) and the b^s* state of expression (18) must interact so as to 
equalize the long-range energies of a ^2+ and ^2“ according to observation (2); hence, 

simple correlation studies automatically impose configuration interactions and also pre- 
dict magnitudes and directions of the corresponding energy shifts. Further details are 
provided in the next section. 

Only one other example is needed to illustrate completely the various considera- 
tions used in determining the 18 o(^p) + o(^P) molecular states and electron configu- 
rations. In contrast with expressions (15) to (17), the configurations 

5 n g : (l,l) a + (0,l) b 



and 

5 n u : (l,l) a + (0,l) b 



contain double the number of subconfigurations necessary to involve every constituent 
atomic determinant. The alternative of two ^n g and two ^II U molecular states is not 
allowed because each would have 50 percent probability of dissociating into normal prod- 
ucts and 50 percent probability of dissociating into ionized products, according to whether 

the electrons in ^Og ^(°u) an( * * or exam P^ e ? are evenly or unevenly split 

between the atoms. Configuration interactions analogous to those of reference 4, in 
which the molecular orbital occupations /a g \2 and (cr u \2 short ranges are replaced 


8 



with the valence bond occupation (cr) 2 at large separations, insure normal products and 
result in the doubling displayed here. 

Another difference between configurations (19) and (20) and the state of 

expression (14) is that the (l,l) a + (0,l)b combination does not select a definite g or 
u symmetry as does (l,l) a + (1,1)^; hence, a 5n state of each symmetry must occur. 
All molecular states with the total magnetic spin quantum number equal to 2 and with 
dissociation into normal atoms have now been derived in agreement with the results of 
group theory (ref. 5, p. 321). 

The 12 remaining molecular electron configurations are obtained in a similar man- 
ner and are listed as follows: 



A 3 S* : 2“l^(o)^ ( Jr+ )^^u^ ir g) + ( 7I u)( ir g)( 7r_ ^ + ( a g)(°u) ( ir+ ) 4 ( ?r *) 2 + (tf + ) 2 (tf“) 4 J* ( 24 ) 
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and 

b 1 ^ : ( a ) 4 ( 77 +) 2 ( 77-) 2 (32) 

Absent from this list are combinations of expressions (5) to (13) which correspond 
to zero magnetic spin quantum numbers on both atoms. Unlike the other spin combina- 
tions, such pairings do not represent pure spin states because of the multiplicity of cou- 
plings and phase relations between the atomic spin vectors as they precess about the 
molecular axis. All the pure singlet states are obtained from combinations of atomic 
states with M s = 1 and M s = -1. 

By way of contrast, the molecular electron configurations used by Vanderslice, 
Mason, and Maisch (ref. 3) are given by 
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* n g : 

(°g) 2 (°u)( 1, u) ( ff g)( n u)( ir g) 

(46) 

x n u : 

( a gj^ (°u) (^u) ( 7r g)( 71 u) 

(47) 


: ^gf^u) 2 ^)^) 2 ^) 

(48) 

: 

( CT g) 2 (4) 2 (^g)(^)fi) 2 

(49) 

» ls g 


(50) 


Configurations (33) to (50) differ substantially from those of the present research 
for two reasons: (1) correlations between the oxygen molecule and its dissociation prod- 

ucts are not included in reference 3 and (2) Vanderslice, Mason, and Maisch identify two- 
electron bonds and antibonds more closely with molecular orbital theory than with valence 
bond theory. Consequently, the two-electron occupation (Og j2 often replaces (cr)2, for 
example, and thus fails to note the important configuration interaction with ^a u )2 that 
distinguishes between neutral and ionized products (ref. 4). The simple configurations of 
expressions (33) to (50) are determined mostly from energy arguments which are seldom 
reliable (ref. 4); they do not imply simpler or more symmetric interaction potentials 
compared with those of expressions (14), (15), and (17) to (32). In addition, the configura- 
tions for a b^Sy, and a in expressions (36), (42), (44), and (48) do 

not yield the correct + or - symmetries. (See appendix A.) 

INTERACTION POTENTIALS 

As explained in reference 4, the energy of each individual subconfiguration in 
expressions (14), (15), and (17) to (32) is computed as if that subconfiguration existed 
alone, the interactions between them serving merely to establish different energy asymp- 
totes for different dissociation products. The subconfiguration energies are weighted by 
the squares of their coefficients and then added to yield the atom-atom interaction poten- 
tials (p . All subpotentials <p^ in each subconfiguration are expressed in terms of the 
exchange energies according to the following relations (ref. 4): 

4>2 ~ ±J2 H ±J (51) 
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03 ~ 01 = Ql ± J 1 ~ ±J 1 


(52) 


and 


<p 4 ~ -2J 2 = -2J (53) 

where the subscripts denote the number of electrons in the bond or antibond being con- 
sidered and Qj is the one-electron Coulomb integral. 

The justification for equation (52) follows from screening arguments and the fact 
that Jj is proportional to overlap integrals between atomic orbitals on different atoms, 
whereas J 2 and Qj are proportional to squares of overlap integrals (ref. 4). Since 
the three-electron configuration described by the molecular orbital determinant 

(a+b a-b | a±b) 

is equivalent to the molecular orbital -valence bond arrangement 
(a b | a±b) 

and since the two valence bond electrons effectively screen the atomic cores from the 
molecular orbital electron, <p 3 should be closely approximated by (pi and only slightly 
perturbed by energies proportional to squares of overlap integrals. At the same time, 
the molecular orbital electron is much less efficient in screening the atomic cores from 
the valence bond electrons; hence, the neglect of their J 2 contributions to <p% is con- 
sistent with the retention of J 2 in two- and four-electron bonds pertaining to smaller 
effective nuclear charges and, therefore, to greater overlap between atomic functions. 
Finally, Qj is neglected in equation (52) because it also relates to overlap integrals 
which are smaller than normal as the result of incomplete screening. 

A typical example of the procedure for finding long-range interaction potentials 
from molecular electron configurations is provided in the following manner by the 5^g 
state of expression (14) and the energies in equations (51) to (53): 

0( 5A g) = ~ Jtt) + - J 7 r - 

+ '|(" Jcr + J l7T " J l7l) + ^(" Jcr " + J lir) 

= -(ja + |j 7 r) (54) 
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Each exchange energy in equation (54) is labeled by the a or y electron occupa- 
tion from which it derives, there being no exchanges between electrons occupying orbitals 
of different geometric orientations because of overlap considerations. Positive signs are 
chosen for two- and three-electron bonds such as (cr) 2 and the first of which 

does not exist in this particular example, and negative signs for antibonds such as 
(°g)( CT U ) and ( ir u)( ir g' ) 2, Since a11 exchange energies are negative by definition (ref. 4), 

the ^Ag state of molecular oxygen is obviously repulsive in its outer regions. 

The remaining 17 interaction potentials are deduced in similar fashion and listed in 
the second column of table I. As required by the discussion following expression (17), 
the energy of the a^2g state is shifted downward from the value 



to coincide with the energy of the ^2 a state given by 

^( 5s u) = - J ° ( 56 > 

For such a shift to be possible, the energy of the configuration interacting b ®2g state 
must lie above that of the a 5 2| state, which it does according to the relation 

^b5 S +) = -2(j 0 + J„) (57) 


and the fact that J a has larger magnitudes than J n because of greater overlap between 

p atomic functions, 
z 

Corresponding to the downward shift of the a ^2+ energy is an upward shift of the 
b^2g energy by an amount determined from the use of equations (55) to (57) in the well- 
known variational principle of quantum mechanics. This principle relates the low-level 
(no configuration interaction) approximations <f> a and (p^ for the molecular states a 
and b to the final energies (denoted by primes) through the secular equation 


0 a " 0 ' 
Hab 


Hab 


= 0 


0b " 0’ 


(58) 


If 4> a , cpfo, and one of the cp' solutions are known, the off-diagonal interaction energies 
H a k can be eliminated in the expansion of equation (58) to yield 
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0b = 0a + 0b - 0a 


(59) 


The substitution into equation (59) of a for state a, b ^Sg for state b, 
and for finally gives 


0’(b 5 S+') = 4>(a 5 £+) + <*>(b 5 S|J - 4>( 5 2^ 


= -2/Ja + |J7r 


(60) 


Similar treatments are required for the triplet and singlet 2 configurations with the 
results 

0’( b 3 S+) = 4 >(a 3 S+) + 4>(b 3 2+) - <^X 3 Sgj 

— “ 2 — 2 J a — 0 

= -2^J(j + Jt7^ (61) 

and 

0'(b ^g) = <p(& X Sg) + 0^b ^ 

= (J a - - 2^J a - - Jqt 

= -2(j CT - | jJ (62) 


Also shown in table I are the Vanderslice, Mason, and Maisch interaction potentials 
(ref. 3) computed from the electron configurations of expressions (33) to (50) and the 
assumption that each electron occupies a molecular orbital. The electron energies are 
defined arbitrarily to be J 2/2 (that is, j/2) in attractive orbitals and -3J2/2 (that is, 
-3J/2) in repulsive orbitals, so that the two-electron subpotentials 


and 



- 2(Jq/2) 


- J a 



(63) 


(64) 
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satisfy valence bond requirements. The four-electron subpotential is given correctly by 
the expression 


0 4 



= 2(J ct / 2) - 2(3J a /2) = -2J a 


(65) 


However, the three-electron subpotentials 


03 


W 


! ( a u)| 


= 2(J a /2) - 3J a /2 = -J a /2 


and 



= Jo/2 - 2 (3 Ja/2) = -5J CT /2 


( 66 ) 


(67) 


calculated by the procedure of reference 3 fail to satisfy equation (52) because no distinc- 
tion is made between exchange energies varying as the first and second powers of the 
overlap integrals. Equations (66) and (67) are nonsymmetric with respect to the dissocia- 
tion asymptote, neither equation predicting an attractive force; hence, stable configura- 
tions such as the helium molecule-ion would not be possible within this framework. A 
further discussion is given in appendix B. 


DISCUSSION OF RESULTS 


Starting with correlations between molecular and atomic states and continuing 
through configuration interactions and the assignments of individual electron energies, 
the methods of the present research and reference 3 are so different that the complete 
lack of correspondence between the long-range interaction potentials listed in columns 2 
and 3 of table I is understandable. Just as in the previous calculations for nitrogen 
(ref. 4), the present modifications greatly enhance the symmetry characteristics of the 
potential energy curves with respect to both the dissociation asymptote and each other. 
Among the useful relations between potentials are the following expressions: 
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( 71 ) 


(72) 


In both sets of calculations in table I the appearance of only two exchange energies, 
J a and J w , permits all 18 long-range interaction potentials to be expressed as linear 
combinations of two given ones; however, even this much complexity may not be required 
for computations of high-temperature transport coefficients. Since overlap considera- 
tions imply magnitudes of J CT far greater than in the range of internuclear separa- 
tions most important for high-temperature transport properties, a good approximation 
should be the setting of all potentials equal to zero except those found from the second 
column of table I to satisfy 

t^AgJ ~ $(& X 2g) « ~ -</}(& 5 2gj ~ -4>( 5 2 u) 

® -<*>( 5 Ag; * -2~ 1 cp(b 5 sj) « -2 -1 ^b 3 rj) ~ -2“ 1 0^b X 2g) « J CT (73) 

The weighted collision cross sections corresponding to these potentials then combine to 
give the effective value 

4ft 94 Q 

Qmean Q( < ^ >a ®) + ~r~~ Q(<£~±Jo) + ~ Q(0 ~-2 J or) 

o 1 81 81 

Q(4>*±Ja) + | Q(<p a -2Jo) (74) 

where A-doubling is included in the coefficients. 


A similar analysis of the third column in table I yields 


26 


Qmean * ff QM^o) 


+ |y Q ~ 2~ 7 + h q(< ^ _2Jct) 


(75) 


which is larger than the effective cross section in equation (74) because of nonzero con- 
tributions from every molecular state. Only 41 percent of the modified potentials con- 
tribute to the effective cross section; hence, transport coefficients based on the results 
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of reference 3 may be somewhat smaller than the true values for a gas of pure atomic 
oxygen. 

Further comparisons between the two sets of long-range interaction potentials and 
collision cross sections are not very productive at this stage because of the lack of good 
transport data at high temperatures. Reliable formulas for J a also are not available. 
However, several interesting conclusions and observations can be made by extrapolating 
and adjusting the molecular electron configurations of expressions (14), (15), and (17) to 
(32) to cover smaller internuclear separations. 

EXTRAPOLATIONS OF SELECTED ELECTRON CONFIGURATIONS 
TO INTERMEDIATE INTERNUCLEAR SEPARATIONS 

Two fundamental difficulties are associated with the preceding analysis: (1) the 
internuclear separation below which the configuration interactions between similar E 
states cease to be significant is not known and (2) the close-range limit of the region in 
which combinations of electron subconfigurations are necessary for correlation purposes 
is uncertain. Whereas the first difficulty is resolved if J n is negligible in equa- 
tions (55) to (62), the second requires considerable thought for most interactions. 

Although some evidence is given in reference 4 that relations between potentials are 
more consistent with intermediate-range calculations and spectroscopic data when corre- 
lation is maintained, other arguments point to the breaking of certain types of correlation 
at relatively far distances. In principle, of course, the statistical weights of the electron 
subconfigurations in any molecular state can be obtained as functions of the separation 
from ab initio quantum calculations, but this approach is beyond the scope of the present 
research. 

One example of the type of correlation which may be broken or altered at fairly 
large separations is displayed by the state of expression (27). A likely evolution 

of the electron configuration as the internuclear separation decreases is given by 

1 A g : 2~ 1 (cr) 2 (it+) 4 (7t-) 2 + (TT+) 2 (n~) 4 + ^uf^gj^uf^g) + (^)(^g) 2 (^u)^g)^ 

— 2" 1 (cr)2 ( 7r+ )4 (7T - ) 2 + (r^V) 4 + 2 ^+) 2 ("gfaf (”gj 

^ 2 " 1 (° r g) 2 f 7r+)4(7T " )2 + (7r+)2 ^" )4 + 2 ( 7r 5) 2 ( 7r i)( ir u) 2 (^)] 

- 2_1 ( ff g) 2 ^ + > 4 ( 7r u) 2 + (0 2 (*-) 4 + 2 ( 7r u) 2 ( 7r g)( 7r u) 2 ( 7r i)J <™> 
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where the second step represents an alteration in the 77-orbital occupations corresponding 
to the lowest energy choice. The reason for this change being likely at relatively large 
distances is the fact that it can occur without suppressing any of the four combinations 
inherent in collisions between (1,1) and (1,-1) atoms; in addition, the configuration inter- 
action between and which is rec l uired at very large 

distances to distinguish between normal and ionized products, should dissipate rapidly 
with decreasing separation. 

Steps 3 and 4 of expression (76) correspond first to the substitution of (°g)^ for 
(ct) 2 and then to the substitution of ( 7 r u ) 2 for ( 7 r) 2 , the priorities relating to overlap 
influences on the dissipation of the two-electron configuration interactions discussed in 
reference 4. Although the limits on the internuclear separation for the various stages of 
expression (76) are not defined by this procedure, the qualitative evolution of the molecu- 
lar electron configuration to the less specific notation ^g^^u^^g^ I s confirmed by 
Herzberg (ref. 5, p. 345). Finally, the energy evolution associated with expression (76) 
is given by 

= Ja - \ Jtt — Ja + Ju - \ J tt 

2J i 0 + - ~ 2J io- + 2JJJJ. - (77) 


a 1 ^ 


The two remaining molecular states evolving to 
and according to the following expressions: 


configurations are 




r 1 / 2 (o ) 2 


( 7 T +) 4 ( 77-) 2 + ( 7 T +) 2 ( 77-) 4 


and 



(78) 


(Equation continued on next page) 
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Corresponding to these configurations are the energy evolutions 



— ^2Jjq- + 2J - 2J 7l (80) 


and 

tf>(x 3 Sg) = 0 — J l7r — J a + 2J U 
2J 1 a + 2J l7 t 


(81) 


where recognition of the interaction with b *2g is taken in the first step of equation (80). 

A convenient summary of equations (77), (80), and (81) near the equilibrium separa- 
tion is contained in the expression 





±4>(x 3 Z^\ + 4 aV 


(82) 


whereby the energy of the *Ag state lies midway between that of X 3 2“ and a *Sg. 

Also established by the same equations are X 2 Zg as the ground state and a^Eg as 

the highest of the three most stable configurations. These results agree perfectly with 
the spectroscopic ordering of equilibrium energies and are surprisingly close to the 
observed ratios (ref. 5, pp. -346 , 446, and 560) 


E a A S 


1t + 


E( i A, 


E(X 3 E" 


= 2.63 


1.98 


1.00 


(83) 


20 



Responsible for this success are the correlation analyses of the preceding para- 
graphs, which trace from the isolated atoms and thus permit a more detailed study of the 
three molecular states characterized by the Herzberg notation Without 

such detail, Herzberg (ref. 5, pp. 358 and 359) concluded that the X 3 S~ energy is 
shifted downward from the other two because of configuration interaction with a 3 2" 

o 

state dissociating into 3 P and 4) atoms. This assumption no longer seems necessary 
to explain the basic energy ordering, although it may well contribute to the asymmetric 
spacing implied by equation (83). 

A second group of three stable configurations corresponds to the Herzberg notation 
(CTg) 2 (7 r u )3(u-g) 3 and to the following evolutions as the internuclear separation decreases: 



and 



— 2 


- 1/2 

- 1/2 


W 2 jl 

[<?( 'gX*; 

+ KX 7 ^) 

( 7r u) 2 ( 7r g) 

N 2 

r +\2 / +\ / 
(%) 

T u)( w i) 2 + (4)(^g) 

2/ _\2 / 
(%) yfg) 


(85) 


( 86 ) 


21 



Associated with these expressions are the energy evolutions 


4>(a 3 2+) = 0 — -2J ff J a - 3J V — 2J la - 3J n (87) 

0(3 A u ) = -J* — J a - | 3 V — 2J la - | J* (88) 


and 

«»( ls u) = J a-2J lff (89) 

so that 



The states A 3 2j, 3 A U , and *2“ are in the order of increasing stability near the 
equilibrium separation of nuclei. 

Since equations (87) to (89) do not involve the contributions to equations (77), 

(80), and (81), the energies of the second group of three molecular states ^i.e. , Iz^, 

3a u , and A 3 2^) lie above those of the first group at equilibrium distances. This obser- 
vation plus the energy ordering of the second group and the relation displayed in equa- 
tion (90) are in general agreement with the conclusions of Moffitt (ref. 6) and Fumi and 
Parr (ref. 7), but not with Herzberg (ref. 5, pp. 446 and 560; also, ref. 8), who lists 3^ u 
as the lowest state in the second group. 

An overall survey of the six stable configurations discussed in the present section 
suggests crossings of the X 3 2" potential energy curve with those of the ^Ag, a *2jf, 

and states and the 3^ u curve with that of A 3 zj; in addition, potential humps or 

maxima located somewhere between equilibrium and large internuclear separations are 
predicted for the 3 a u and A 3 2* states. Potential humps also may occur in the ^Ilg 

and 3 n u curves, all n states being at least slightly stable at sufficiently small 
separations. 

3 _ 

Another interesting observation concerns the fact that the ground X Zg state is 
not the most attractive one at large distances, that distinction falling to a ^Z| and ^Z^; 

meanwhile, the most repulsive configuration belongs to b 3 Zg. Few of these conclusions 
are shared by Vanderslice, Mason, and Maisch (ref. 3), who also fail to distinguish 
22 



between the individual energies in each triad of strongly stable molecular states. Many 
important details are either missing or incorrectly described in reference 3 because of 
the lack of correlation between molecular and atomic configurations and the arbitrary 
assignments of molecular orbital energies. 

Finally, the energy evolutions of equations (77) and (80) forecast a much better suc- 
cess for the extrapolations of Rydberg-Klein-Rees semiempirical *Ag and a *2* 

potentials to large distances than for the extrapolation of the curve in equa- 

tion (81). As observed in reference 3, the outer extremity of the Rydberg-Klein-Rees 

3 - 

X Sg curve is too steep to be connected smoothly with Vanderslice, Mason, and Maisch’s 
strongly attractive valence bond energy at still larger separations. A rather sharp 
bending of the Rydberg-Klein-Rees curve leading to a very gentle approach to the dissoci- 
ation asymptote is implied; hence, a good indirect confirmation of the flatness of the long- 

Q _ 

range X J Sg tail in the present analysis is obtained. 

CONCLUDING REMARKS 

Correlations between the electron structures of oxygen molecules and their dissoci- 
ation products are shown to significantly alter previous determinations of the 18 long- 
range interatomic forces operating in a pure gas of normal oxygen atoms; in addition, 
certain correlations predict and specify the energy shifts resulting from configuration 
interactions with states of the same symmetry. Unlike much of the existing body of 
research, a consistent mixture of valence bond and molecular orbital energies is assigned 
to the various bonding and antibonding configurations. Reasonable extrapolations of the 
electron structures and potential energies to intermediate distances successfully predict 
the ordering of the six most stable molecular states and suggest a number of potential 
humps and curve crossings in the process. 

Among the conclusions pertaining to large internuclear separations is the predic- 
tion that the transport coefficients for atomic oxygen have been underestimated at high 
temperatures Uy researchers using the approximation of perfect pairing. Because of the 
correlation analyses and other considerations in the present paper, nine of the 18 atom- 
atom interaction potentials are nearly flat in the outer regions and, thus, do not contribute 
effectively to the collision integrals and cross sections of kinetic theory. Some confirma- 
tion of this flatness is obtained for the ground state of molecular oxygen by studying the 
long-range extrapolation of the corresponding Rydberg-Klein-Rees semiempirical poten- 
tial energy curve. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 5, 1970. 
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APPENDIX A 


ELECTRON CONFIGURATIONS AND INTERACTIONS 

To a large degree, both the present research and that of reference 3 involve the 
writing of valence bond energies directly from simple molecular orbital descriptions of 
electron configurations. However, in expressing double cr-orbital occupations by the 
symbol (cr) 2 instead of or ( a u)^ f° r electrons with paired spins, for example, 

the present paper explicitly recognizes a major difficulty with molecular orbital language 
that reference 3 does not — namely, that configuration interaction between and 

^a u ^ 2 at large internuclear separations makes the atomic orbital or valence bond descrip- 
tion more appropriate. The statements in reference 3 that a pair of electrons in a bonding 
molecular orbital contribute the valence bond energy J and that the aforementioned con- 
figuration interaction is a second-order effect are individually incorrect and mutually 
inconsistent. Molecular orbitals cease to exist when such configuration interactions are 
important, so that the practice of arbitrarily assigning valence bond energies to these 
configurations cannot possibly give relations between the two descriptions. Extensions of 
this procedure to three-electron bonds and antibonds are certain to yield significant 
errors. (See the text and appendix B.) 

Another difficulty with the notation and electronic states of reference 3 is the appar- 
ent confusion in some expressions of the meaning of the + and - superscripts on the 
77 -orbitals. Vanderslice, Mason, and Maisch clearly adopt the convention used in the 
present paper of letting these superscripts represent the reflection symmetries of the 
orbitals with respect to a plane containing the molecular axis, which is by far the most 
convenient interpretation for the geometrical and overlap considerations so essential for 
the subsequent applications of the approximation of perfect pairing. In some states, how- 
ever, these authors seem to inadvertently change this convention to one in which the same 
superscripts refer to the two different senses of rotation of a single electron about the 
molecular axis. For example, their configuration for the a ls+ state can be expressed 
( 77 -electrons only) as 

^u + u" g + |u + u" g") 

which does not produce a 2 + state unless the orbitals are defined by the angular momen- 
tum eigenfunction relations 

U± = 2" 1//2 ^7T+ ± iir^j (Al) 
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and 


g± = 2 


- 1/2 



(A2) 


The left-hand and right-hand sides of equations (Al) and (A2) invoke the two different 
superscript conventions. 

According to references 6 and 7, the complete wave function for the a state 
near the equilibrium separation is given by 

2” 1 /‘ i [ju + u- g + |u+ u- g") + (u+ u- g~|u + u- g+) 
so that the electron configuration becomes 




(A3) 


with the aid of equations (Al) and (A2). A second •'■Eg state considered in these refer- 
ences is 


lv+ 

S g 


,- 1/2 


(u + g+ g- | U" g + g“) + (u- g + g" I u + g+ g") 


2 " 1/2 (°gy 


(77+) 4 (77g) 2 + 


( 7 r gf ( 7r ") 4 


(A4) 


the interaction of which with expression (A3) gives 

a J Sg : 2“ 1//2 ^CT g ^ 2 |(7r + ) 4 (T7") 2 + (7r l ’) 2 (Tr")^J (A5) 

for the inner regions of the potential tail. 

Hence, the configuration interaction which Vanderslice, Mason, and Maisch neglect, 
at least verbally, is responsible for the destruction of 77 -molecular orbitals and the intro- 
duction of valence bond concepts in a completely systematic way. The arbitrary assign- 
ment of valence bond energies to molecular orbitals can only lead to confusion and incor- 
rect extrapolations to more complicated systems; it can never explain the basic process. 

Expressions (A3) and (A5) are identical with the third and second steps, respec- 
tively, of expression (78), the first step of which corresponds to an additional configuration 
interaction to replace ( a g) 2 with (a)2 at larger internuclear separations. Molecular 

orbitals persist at longer ranges for a-electrons because of greater overlap. The vast 
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differences between the electron configuration of expression (A5) and that given in refer- 
ence 3 for the a ^Eg state are caused primarily by the aforementioned misinterpreta- 
tion of + and - superscripts on 7r-orbitals. 

Similar studies of the bound 1 Ag, X 3 E”, A 3 Su, 3 A U , and *E{J states yield the 

same conclusions and show configurations based on references 6 and 7 to be identical with 
one or more steps of expressions (76), (79), (84), (85), and (86). 
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APPENDIX B 

THREE -ELECTRON ENERGIES 

As mentioned in the text and discussed more fully in reference 4, Vanderslice, 
Mason, and Maisch (ref. 3) have proposed a formal technique for extrapolating two- 
electron energies to three-electron configurations. The essential feature of the method 
is the assignment of valence bond energies to single electrons, all reference to electron 
pairs being eliminated. Although correct energies can be obtained in this manner for 
two- and four-electron bonds and antibonds, the procedure is artificial because it over- 
looks in its verbal description, but not in its energy assignments, the important correla- 
tion effects that keep electrons apart. Hence, the same energies applied to singly 
occupied molecular orbitals (as in one- or three-electron a and ir configurations) 
impose correlation effects which do not exist in such examples. 

The mathematical and conceptual consequences of imposing nonexistent correlation 
effects on the electron energies of singly occupied molecular orbitals are severe. No 
longer can the helium molecule-ion form a stable state, nor do the long-range attractive 
and repulsive potential tails of the hydrogen molecule-ion appear as approximate mirror 
images of each other across the dissociation asymptote if the interaction energies J/2 
and -3J/2 of reference 3 are assigned to individual electrons. Magnitudes of molecular 
orbital interaction energies are changed from dependencies on overlap integrals between 
atomic functions to dependencies on squares of overlap integrals by the introduction of 
nonexistent correlation effects. Theories which ignore such distinctions cannot possibly 
apply to molecules (for example, oxygen) in which not all a and v configurations are 
comprised of even numbers of electrons. 

True correlation effects in two-electron bonds are gradually dissipated as the inter- 
nuclear separation is decreased from the regions of the potential tails considered both in 
this paper and in reference 3. Simultaneously, the configuration interactions that cause 
the correlation effects and produce (o)2 from (°g)^ and f° r example, also dis- 

sipate and thereby lead to the full recovery of molecular orbital descriptions. Each elec- 
tron now traverses the molecule more or less independently of the other and, except for 
the comparatively small Coulomb interactions between them, can be assigned single- 
electron molecular orbital energies. Only for such intermediate-to-small internuclear 
separations is the one-electron bond strength approximately half that of the two-electron 
system. Vanderslice, Mason, and Maisch’s extrapolation of this result to larger separa- 
tions is not valid because the advent of configuration interactions and associated correla- 
tion effects makes the two-electron bond energies proportional to squares of overlap 
integrals instead of the direct proportionalities characteristic of one- and three-electron 
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configurations; hence, the two-electron bond is actually weaker than the one-electron bond 
in the regions of potential tails. Crossings of one- and two-electron bond energies must 
occur as the distance increases from equilibrium separations. 

Two-electron antibonds, on the other hand, are always best described by valence 
bond theory because of the equivalence of 

(a+b a-b) 


to 


(a b) 

Correlation effects persist in this example at all internuclear separations; consequently, 
the assignments in reference 3 of different energies to the electrons presumably occupying 
different molecular orbitals is fictitious. Neither electron acts independently of the other 
and therefore the behavior and properties of neither electron can be compared with the 
situation prevailing when the other electron is absent. 
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